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ABSTRACT 

We report on the statistics of strong (Wr > 0.15 A) C IV absorbers at 
z = 1.5 — 3.5 toward high-redshift gamma-ray bursts (GRBs). In contrast 
with a recent survey for strong Mg II absorption systems at z < 2, we find 
that the number of C IV absorbers per unit redshift dN/dz does not show a 
significant deviation from previous surveys using quasi-stellar objects (QSOs) 
as background sources. We find that the number density of C IV toward 
GRBs is dN/dz\G^B{z ~ 1-5) = 2.2+^-8, dN/dz\GKB{z ~ 2.5) = 2.3+};^ and 
dN/dz\Q^-Q{z ~ 3.5) = l.l;!;Q;g. These numbers are consistent with previous 
C IV surveys using QSO spectra. Binning the entire dataset, we set a 95% c.l. 
upper limit to the excess of C IV absorbers along GRB sightlines at twice the 
incidence observed along QSO sightlines. Furthermore, the distribution of equiv- 
alent widths of the GRB and QSO samples are consistent with being drawn from 
the same parent population. Although the results for Mg II and C IV absorbers 
along GRB sightlines appear to contradict one another, we note that the surveys 
are nearly disjoint: the C IV survey corresponds to higher redshift and more 
highly ionized gas than the Mg II survey. Nevertheless, analysis on larger statis- 
tical samples may constrain properties of the galaxies hosting these metals (e.g. 
mass, dust content) and/or the coherence-length of the gas giving rise to the 
metal-line absorption. 
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Introduction 



After extensive studies on absorption line systems in spectra of quasi-stellar objects 
(QSOs) it has been widely accepted that the majority of systems discovered are intervening 
and that they trace the cosmological expansion of the universe. Various surveys have been 
conducted that focus on the cosmological evolution of metal absorption systems for distinct 
atomic transitions. A standard observational measure is the number of absorbers per unit 
redshift, dN/dz, which includes the cosmological evo lution in the physical pathlength and 



evolution intrinsic to the absorbers themselves (e.g., iLanzetta et al. 1119871 : ISargent et al. 



1988l : ISteidellll990l ). 



Since the production of heavy elements takes place in the stars which mainly group 
into galaxies, it is natural to relate metal ab sorbers to galaxies (although this has prove n 



difficult to test, especially at high redshift; e.g.. I Churchill et al. 112004 : iTripp &: Bowen 



20051 ). 



Therefore, studying the evolution of metal absorbers may impact our understanding of both 
galaxy evolution and the physics of the intergalactic medium (IGM). 



Recently, iProchter et al. I (j2006bl ) have found a significant overabundance (a factor of ~ 
4) of strong (rest-frame equivalent width of Mg II 2796, VFr(2796) > 1.0 A) Mg II absorbers in 
lines-of-sight toward gamma-ray bursts (GRBs) when compared to the statistics drawn from 
lines-of-sight toward QSOs. This result is striking because a key hypothesis of this experiment 
is that intervening absorbers are independent of the background source. Several physical 
effects have been proposed to explain the overabundance: (1) dust within the Mg II absorbers 
may obscure faint QSOs, (2) the Mg II gas may be intrinsic to the GRBs, (3) the GRBs may 
be gravitationally lensed by these absorbers, and (4) the absorbers are s mall enough that 



different 'beam sizes' between GRBs and Q S Os may affect the statis tics (IProchter et al 



2006bl : iFrank et al. I l2006l : iHao et a/. 1120071 ). IProchter et al I (l2006br ) and IPorciani et al 



(120071 ) have argued that none of these explanations is likely to explain the full effect but it 
may be possible for several to contribute together to resolve the discrepancy. 

In this paper, we revisit high-resolution GRBs spectra and look for C IV absorbers to 
obtain the first statistics of such absorbers in this type of lines-of-sight. Because the C^'^ 
ion has a much higher ionization potential than Mg"*", the C IV doublet is likely to trace 
more diffuse and hotter gas. Therefore, C IV systems may represent a different "population" 
of larger "cross section " absorbers and the stu dy of their statistics may help address the 
problem opened by the IProchter et al. I (l2006bl ) result. 



The paper is organized as follows. In § [2] we present the spectra used in this study. 
The statistical sample and the redshift number density calculation are described in § |3l The 
results are summarized in § HI and discussed in § O 
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Data 



We have drawn our spectral sample from echelle observations summarized in Table [TJ 
We have restricted our analysis to the sightlines with at least coverage redward of the Lya 
forest (i.e. Zq'^q < z^^^ < Zq^j^^) and to spectra with moderate to high signal-to-noise ratio 
(> 5pix^^). The sample is comprised of echelle spectra of the following GRB afterglows: 
GRB021004, GRB050730, GRB050820 GRB050922C and GR B060607. The data was ob - 
taine d with the MIKE/L as Campanas jBernstein et al. I l2003h. UVES/VLT JPekker et al. 



2 pool ) and HIRES/Keck (IVogt et al. Ill994j ) echelle spectrographs allowing high spectral res- 
olution of -R ~ 40000. The UVES data were retrieved from the ESQ archive and processed 
with the standard UVES pipeline. For each GRB the individual spectra acquired with a 
given instrument were reduced, coadded and normalized with standard te chniques. We refer 



the reader to the following papers for a ful l discussion of the observations (IFiore et al. 



Chen et al. 2005: Prochaska et al. 2007: Piranomonte et al. 2007; Ledoux et al. 2006). 



2005 



3. C IV Absorption Systems toward GRBs 

3.1. Equivalent Widths 

Equivalent widths were calculated using both direct pixel integration, W'^^^'^\ and a 
Gaussian fit, 14^5"^"*'*, applied directly to the high-resolution data (since previous surveys 
used spectra of lower resolution, we first confirmed that smoothing and re-binning our data 
does not alter the equivalent widths). We obtained VI^p*^'^' over the smallest velocity window 
that contain the whole system. The rest frame equivalent width Wr was calculated by taking 
Wr = Wobs/i^ + z). We arbitrarily defined the redshift of the absorption system such that 
it appears centered in velocity space. 

To obtain W^""^^^, we fitted a single or multiple Gaussian profile to the absorption lines, 
and summed each component of an individual system. In the case of multiples components, 
the error associated is, cr'^ystem — Xli '^i^ where cxj is the individual error of component i. 
Figure [1] shows that equivalent widths from both methods are, in general, self-consistent. 
The outlier is due to the blended feature in the main component of the system at z = 2.89048 
(see Figure [2]). Hereafter, we will use only W^"'^^'^ because of its robustness (continuum and 
blended systems are fitted, independently of velocity window). 
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3.2. Definition of the Samples 

3.2.1. Full Sample 

We first define a Full Sample of C IV absorption systems that satisfy the following 
criteria: (1) ^^iJB < ^Ss^ — ^grb'i i-^-' Full Sample includes lines associated with the 
GRB host galaxy (we will later restrict the analysis to intervening C IV systems). (2) The 
equivalent width of the C IV doublet must be detected at the 5ct level or higher. No conditions 
were imposed to the equivalent width ratio for both components although all members of the 
Statistical Sample (see below) exhibit an equivalent width ratio iy(1548)/iy(1550) between 
1 (saturated regime) and 2 (unsaturated regime). (3) A complex system is considered a single 
system if the velocity components group within 500 km s^^. This last condition is useful 
to compare w ith previous surveys that were made at lower spectral resolution (for instance, 



Steidel II1990I . used a velocity window of 150 km s^^ to group components of a single system 
but then defined a statistical sample counting as a single system all components within 1000 
km s~^; Our Statistical Sample, defined below, does not have systems separated by less than 
1000 km s~^, therefore these two samples are comparable). 

Although performed redward of the Lya forest, the C IV survey has one significant 
contaminant: the O I 1302/Si II 1304 pair of transitions has nearly identical separation as 
the C IV doublet. It is trivial, however, to identify these contaminants by searching for other 
Si II transitions and the C II 1334 transition. 

We found 29 candidate C IV absorption systems in our Full Sample. This number 
includes absorption systems with z ~ zgrb- Table [2] shows the information on redshift 
limits of each line-of-sight, z'^^^ and Wr- 



3.2.2. Statistical Sample 



To compare the statistics of C IV absorption along GRB sightlines with that one for QSO 
sightlines, one must define a sample with identical equivalent width liraits. Previous surveys 
have defined this limit to Wmin = 0.15 A in the rest frame ( Steidel Ill990l : iMisawa et al 



2OO2I ) . We adopt the same value, i.e., Wy. > 0.15 A in both members of the C IV doublet. 



In addition, to avoid contributions from galaxies clustered w ith the hosts of G RB and 
wind features that might be associated with the GRB progenitor (IGhen et al. 1120071 ) and to 
compare with previous surveys, we define a Statistical Sample with t he same Zstart as in the 
Full Sample but Zend at 5000 km from zgrb (as in lSteidel Ill990l . for zqso)- 
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All the spectra included in our Statistical Sample permit detections of lines with Wr > 
0.15 A with at least 5a significance from to 5000 km within Zq^j^. The Statistical 
Sample is composed of 7 absorption systems (see systems labeled with a dagger in Table [2] 
and see Figure [2] for velocity plots). We excluded 2 systems due to their dubious character 
(systems labeled with an 'X' in Table [2]): Zabs = 3.254, although not saturated, it shows an 
equivalent width ratio < 1; Zabs = 2.3589 falls very close to Wmin A but W/''"""(1550) < 0.15 
A. 



3.3. Redshift Number Density 

The redshift number density, dN/dz, is the number of absorbers per unit redshift. By 
definition, dN/dz = Nabs/ ^z where Nabs is the number of absorption systems in the redshift 
path Az. Defining the number of lines-of-sight between z and z +dz as g{z) (a function that 
defines the redshift path density), the redshift path can be written as = Yli9{.z)^z where 
5z is the grid resolution. Figure [3] shows g{z) for our sample. Figure H] shows our result on 
dN/dz toward GRBs {dN / dz\QYiB) ■ The choice of redshift bins is arbitrary; we used 3 bins 
of 1 redshift unit between z = 1 and z = 4. 

To estimate the error in dN/dz we assume a Poissonian distribution in the number 
of absorbers Nabs in a given redshift bin Az**" and then assign the upper and lower limit 
corresponding to la confidence le vel of a Gaussi an distribution, (Icr = 0.8413) lo"iv^^, and 

dN/dz ^ ^'^NatJ 



Vat, 

is done for each redshift bin. 



Icr^ respectively obtained from iGehrels I (119861 ). Therefore, cr^^,^^ = la^' /A^;^*". This 



4. Result 



For the Statistical Sample we find dN / dzlc^^iz ~ 1.5) = 2.2^^;4, dN / dzlc^B^z ~ 2.5) 



2.31};^ and dN/dz\GRB{z 
bin. 



3.5) = l.l_o;9. Table[3]shows the specific results for each redshift 



How do these numbers compare with dN/dz\o^ci7 One of the l argest published sur- 
veys of C IV along QSO sighthnes was carried out by ISteidel I (Il990l ). In that survey, 66 
QSO spectr a with low spe c tral r esolution (~ 1 — 2 A) were used (55 of them w e re ob 



tained from Sargent et al. (jl988 )). We used their result as comparison. Steidel (1990) 



measured the incidence of C IV systems with W„ 
as dN/dz\Qso ■ 



0.15 A and parameterized their results 



Steidel 


(1990) and 


Sargent et al. 


IIUVVD, 

(1988) 
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published survey by iMisawa et al. I ( 20021 ) (18 QSO spectra with resolution of ~ 2 A) al- 
though this survey has much lower redshift coverage and does not provide a better sensitivity 

(iy^,„ = o.i5 A). 

We find that dN / dz\Q-Rs, decreases with redshift in the range between z = 1 and z = 4 
as in QSO surveys. In Figure H] we compare dN/dz\GKB with dN/d z\qfi,o using our binning. 
It can be seen that, within errors, our result on C IV matches the ISteidel I (Il990l ) and the 
Misawa et al. I (120021 ) ones (according to Misawa et al, the fact that these two previous 
results do not match each other would be due to an statistical accident). On the other 
hand, if we use only one redshift bin between z = 1 and 2; = 4 we find that the incidence of 
C IV absorbers along GRB sightlines is less than twice that one along QSOs (Steidel 1990 
sample) at the 95% confidence level, i.e., a large overabundance is improbable (even taking 
into account the 2 dubious systems excluded from the Statistical Sample). Therefore, there 
is no significant difference between the GRB and QSO statistics in our sample. 



Discussion 



The principal result of our analysis is that the incidence of C IV absorbers {Wr > O.ISA) 
along GRB sightlines is consis tent with that observed a long QSO sightlines. Note that the 
same conclusion is reached by ISudilovsky et al. I (120071 ) using similar data (although these 
authors use a smaller total redshift path and a lower equivalent width cutoff than ours). 

Quantitatively, we set a 95% upper hmit of dN/dz\^^^ < 2dN/dz\QgQ over the redshift 
interval 2; ~ 2 to 4. At the surface , this result lies in stark contrast to the incidence of strong 
Mg II absorbers at 2; < 2 where IProchter et al. I ( 2006bl ) find a lower limit dN/dz\^^^ > 
2dN/dz\^^Q at the 99.9% significance level. 

Are these two results contradictory? Let us first stress that the two absorption-line 
samples may be nearly disjoint. Firstly, one important difference is that the Mg II absorbers 
were surveyed at z < 2 while most of the pathlength surveyed here for C IV absorption 
has z > 2. While there is no indication that the dN/dz\^^^ enhancement is declining with 
increasing redshift, there are possible reasons to expect such an effect (e.g. gravitational 
lensing efficiency is maximal at ziens ~ zgrb/'^)- In that case the Mg II enhancement would 
vanish at the redshifts probed here with C IV. 

But even if this were not the case, let us consider a second argument in favor of the 
present results. The key issue is that C IV may track a distinct phase of gas from Mg II in the 
intergalactic medium owing to its substantially higher ionizatio n potential. This is refl ected 
in the statistics of QSO surveys. At 2 = 2, dN/dz\^fQ ^ 0.4 (IProchter et al. Il2006ah and 
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dN/dz\^^o ^ 2.5 flSteidellll990[ ). i.e., C IV systems factor of 6 more abundant than 

Mg II for the Wr thresholds relevant to our study. In other words, there is plenty of C IV 
(weak) systems that do not show strong Mg II absorption. 

On the other hand, we expect that nearly every strong Mg II absorber will also exhibit 
C IV absor ption. Furtherrnore, m any of these will have an equivalent width in excess of 
O.ISA (e.g. Churchill et al. 1999 ). Therefore, if the enhancement in dN/dz^^^ continues 
beyond 2; = 2, we do expect a certain bias to larger dN/dz\^^-Q. However, since there 
is not a one-to-one correspondence between Wr > lA Mg II absorbers and Wr > O.ISA 
C IV absorbers, such possible bias is diluted in the much more numerous statistics of C IV 
(indeed, none of the C IV absorbers in our sample that have coverage of the Mg II doublet 
shows significant absorption; e.g., zciv = 1.568,1.989 from GRB 050922C). Even if we 
assume that dN/dz^^^ = 4 dN/dz\Q^Q and that each system also has Wciv > 0.15A the 
effect on dN/dz\'^^^j^ is less than 50% (at most {dN/dz\'^%)/{dN/dz\^^Q) = 9/6 given that 
{dN/dz\S)/idN/dz\^i^) = 4/1 and {dN/dz\'^^^)/{dN/dz\^l^) = 6/1). 

While the change to dN / dz\^^-Q for ^^1548 > 0.15A may be modest, as argued above, it is 
possible that an enhancement in dN/dz\^^^ will imply an C IV equivalent width distribution 
that is very different from that toward QSOs, e.g., a much higher incidence of pyr(154 8) > lA. 



Figure [5] presents a histogram of the equivalent widths for our analysis and that of ISteidel 



(I1990I ) with arbitrary binning. We performed a K-S test on the un-binned distributions of 
rest frame equivalent widths that rules out the null hypothesis (of similar distributions) at 
only the 1% c.l. Therefore, there is no indication of possible differences between equivalent 
width distributions of C IV absorbers toward GRBs and QSOs. 

In summary, there is no fundamental conflict between observing an enhanced incidence 
of strong Mg II absorbers along GRB sightlines without a corresponding enhancement of 
C IV absorbers. Before concluding, let us speculate on the implications of these results 
under the assumption that future data confirm the currently measured incidences. If our 
results on dN/dz\'Q^^ and the current enhancement in dN/dz\^^^ are confirmed in larger 
surveys, one may use these contrasting results to infer characteristics of the C IV absorbers. 
For example, if the dN/dz\^^^ enhancement is to be explained by the geometric differences 
between the GRB a nd QSO beams (comparable with the characteristic size of Mg II 'clouds') 



(IFrank et al. II2006I ). then one would conclude C IV absorbers are significantly larger than 



Mg II clouds and the QSO beam size. Indeed, such differences have been inferred from 



comparisons of metal-line absorption along gravitationally lensed QSOs (e.g.. iLopez et al. 



19991). 



Similarly, the results could constrain the dust redenning or mass of the galaxies hosting 
C IV absorbers. Correlations between these quantities and Mg II equivalent width have been 



- 8 - 



found by rece nt studies using Mg II statistics derived from Sloan Digital Sky Survey QSOs. 



For instance, iMenard et al. I ( 120071 ) h ave parameteri z ed th e (significant) amount of QSO 



redenning due to Mg II absorbers, and iBouche et al. I (120061 ) have found an anti-correlation 



between Mg II equivalent width and mass of the dark-matter halo hosting the absorbing 
galaxy. The different behavior of C IV and Mg II, as revealed by the statistics toward GRBs, 
then can put constraints on these absorber properties by studying Mg II samples selected by 
CIV. 
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GRB 



^GRB 



^start 



021004 



050730 



050820 



2.335 



3.97 



1.715" 



2.965 



2.6147 1.833 



050922C 2.199 



1.511 



060607 



3.082 2.207*" 



2.279" 



3.887 



2.554 



2.146 



3.014'' 



Table 2: Pull Sample: All C IV absorption systems found between -2^^^ < z'-^^^ < z^j^^ for 
each GRB spectra. " The region 1.988 < z < 2.014 is excluded due to insufficient signal to 
noise ratio. The region 2.616 < z < 2.668 is excluded due to a gap in the spectrum, f 
Absorbers in the Statistical Sample, x Systems that were excluded (see section 3.2.2). 
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Redshift Bin 


Nabs 


Az 


dN/dz\^^^ 




[1,2[ 


2 


0.93 


9 9+2.8 


9 7+0.5 
^- ' -0.4 


[2,3[ 


4 


1.74 


2.31}-^ 


9 c;+0.4 
^•'-'-0.4 


[3,4[ 


1 


0.90 


1 1+2-6 

-"-•-■--o.g 


1 1+°-^ 
-'-•-L_o.4 


[1,4[ 


7 


3.57 




9 4+0.3 
^•^-0.2 



Table 3: Specifications of tlie resulting dN/dz\'Q^^ for each bin from the Statistical Sam- 
ple. As c ompa rison the fifth column shows the resulting dN/ dz \ q^q for same binnings from 
Steidel I f ll99oh results. 
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Fig. 1. — Rest frame equivalent width of C IV in ttie Full Sample measured with Gaussian 
fit and with pixel integration. The dashed line is the identity function. 
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Fig. 2. — Velocity profiles of C IV absorption systems toward GRBs in the Statistical Sample. 
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Fig. 3. — Number of GRB lines-of-sight available for a C IV survey with Wr > 0.15A using 
/9c — 5000 km s~^ to discount the contamination of the source. 
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Fig. 4. — Distribution of the number of C IV absorption systems per unit redshift in our 
GRB Statistic al Sar nple. As a cornparison (dashed hne) we show results from QSO survey by 



Steidel I ( 1l990l ) and lMisawa et al. I (120021 ) using the same binning as that for the GRB-C IV 



analysis. The bins have been slightly offset in redshift for clarity. 
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Fig. 5. — Distributi on of the rest fr ame equivalent widths 1^^(1548) for our Statistical Sample 
compared with the ISteidel I (119901 ) one (dashed line). The bins have been slightly offset in 
Wr{lMS) for clarity. 



